Many-body quantum quench in an atomic one-dimensional Ising chain 
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We study non-equilibrium dynamics for an ensemble of one-dimensional atomic Bose-Hubbard chains after 
a sudden quench to the vicinity of the transition point of the Ising paramagnetic to anti-ferromagnetic quantum 
phase transition. The quench results in coherent oscillations for the orientation of Ising spins, detected via 
oscillations in the number of doubly-occupied lattice sites. We characterize the quench by varying the system 
parameters. We find clear evidence for few-body and many-body effects in the oscillatory response. 
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Ultracold atomic ensembles confined in optical lattice po- 
tentials have proven to offer unique access to the study of 
strongly correlated quantum phases of matter |lj i2j- Un- 
precedented control over system parameters as well as excep- 
tionally good isolation from the environment allow for im- 
plementation and quantitative simulation of lattice Hamilto- 
nians |3, 4J, not only bridging the fields of atomic and con- 
densed matter physics in the study of ground-state phases, 
but also opening fundamentally new opportunities to explore 
out-of-equilibrium physics in essentially closed quantum sys- 
tems |l5]|6l. For example, the rapid time-dependent control 
available over system parameters makes it possible to observe 
dynamics arising from a quantum quench, where a parame- 
ter such as the lattice depth is changed suddenly in time 13- 
|9]. Recently it was demonstrated that ID chains of bosonic 
atoms with a superimposed linear gradient potential exhibit a 
quantum phase transition to a density-wave-ordered state, in 
which empty sites alternate with doubly-occupied sites ("dou- 
blons"). Beginning in a Mott Insulator phase of the Bose- 
Hubbard (BH) model ifTOlfTTj . where the on-site interactions 
dominate over tunneling and the atoms are exponentially lo- 
calized on individual lattice sites, a gradient potential is added 
until the potential difference between adjacent sites matches 
the on-site interaction energy, and atoms can again resonantly 
tunnel. This was monitored for individual ID chains with a 
length of about 10 sites with the quantum gas microscopy 
technique Iil2i . and effectively maps onto a ID Ising model 
ifTSl . making it possible to simulate the transition from ID 
paramagnetic (PM) spin chains to anti-ferromagnetic (AFM) 
spin chains in the context of ultracold atoms. 

In this letter, we explore the dynamics of a quantum quench 
for bosonic atoms in such a tilted optical lattice lfT4l [Tsl . 
Specifically, we quench the strength of the tilt to be near 
the phase transition point between PM and AFM regimes, 
and hence take the system far out of equilibrium, inducing 
strong oscillations in the number of doublons, which we de- 
tect through molecule formation. We find clear indications for 
the many -body character of the ensuing dynamics. 

We consider an atomic ensemble in ID in a tilted optical 
lattice potential near zero temperature. For sufficiently weak 
interaction energy the system is described by the ID single- 
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FIG. 1: (color online) Adiabatic tilt through the phase transition 
point. The number of doublons Wj (circles) as detected through Fesh- 
bach molecule formation and the number of singly occupied sites A's 
(triangles) are plotted versus the tilt Ef at the end of the adiabatic 
ramp for V^ = IOEr and a^ = 225(5) ao, giving U = 0.963(20) kHz 
for Vy^j. = 20 Er. The full ramp takes 40 ms. The solid lines are fits 
to guide the eye using the error function. The dashed line indicates 
£c = 1.010(20) kHz. The insets show a schematic of the Ising phase 
transition for a ID Mott-insulator chain from the PM (left) to the 
AFM (right) ground state {lT\. 



band BH Hamiltonian |[TOl augmented by a tilt lfT3l[T6l 

H = jY, a]aj + Y^—hi {hi - \) + EY^ini +Y,eihi . (1) 

('■./> '■ '■ '■ 

As usual, flj (di) are the bosonic creation (annihilation) oper- 
ators at the /th lattice site, n, = djdj are the number operators, 
/ is the tunnel matrix element, and U is the on-site interaction 
energy. The linear energy shift from site to site is denoted by 
E, and e, accounts for a weak external confinement. For suffi- 
ciently small tilt (£ <C U) and sufficiently strong interactions 
(L'^ ^ J) the ground state of the system for one-atom commen- 
surate filling is a Mott insulator with exactly one atom per site. 
When the tilt is ramped from E -^U to E ^ U with E k,U 
at finite 7, the system establishes a regular periodic pattern of 
dipole states, for which one finds double occupancy at every 
second site with empty sites in between llT3l[T4l . see insets to 
Fig.fT] The formation of such density-wave ordering is due to 
an effective nearest-neighbor constraint, which results in cor- 
related tunneling of all particles in the ID chain. When, in 
the course of the ramp, the tilt compensates for the interaction 
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FIG. 2: (color online) (a) Reversibility of the phase transition. The 
number of doublons N^ is plotted as the system is driven from 
the PM into the AFM ground state (f < 40 ms) and back {t > 
40ms) for V, = IOEr. The dashed line is a fit using two con- 
catenated error functions to guide the eye. For this measurement 
Vx.y = 25 Er. (b) Phase transition for various values of the lat- 
tice depth Vj = 8 (circles), 10 (triangles), 12 (squares), and 14Er 
(diamonds). The ramp rate is £ = (17.3, 17.3,8.6,5.8) Hz/ms for 
Vz = (8, 10, 12, 14) Er. The solid lines are fits using the error func- 
tion to determine the half width w. For this measurement Vx.y = 
30 Er. We determine (/ to (1.233, 1.324, 1.401, 1.467) kHz and £c 
to (1.309,1.371,1.431,1.487) kHz for V, = (8, 10, 12, 14)Er. The 
positions of £c are indicated by the dashed lines. The different data 
sets have been recorded on different days, and the variation in N^ 
reflects the day-to-day variations in total particle number (c) Half 
width w of the phase transition obtained from the fit as a function of 
v.. The dashed line shows the enhanced oscillation rate 47. 



energy {E k. U), the particles are allowed to resonantly tunnel 
onto the neighboring sites, however, only when the neighbor- 
ing particle has not yet tunnelled itself. This system can be 
mapped onto an effective Ising spin model lfT2l[T3l . where the 
two distinct ground states (PM and AFM order, respectively) 
are connected via a quantum phase transition with the quan- 
tum critical point £c = t/ + 1 -857 1 17 1. 

To prepare the ensemble of ID Bose chains we start with a 
3D Bose-Einstein condensate (BEC) without detectable non- 
condensed fraction of typically 8.5 x 10* Cs atoms in the en- 
ergetically lowest hyperfine ground state |F = 3,m/r =3) con- 
fined to a crossed dipole trap. Trapping and cooling proce- 
dures are described in Refs. ifTSl [T9JI . The sample is levi- 
tated against gravity by a vertical magnetic field gradient of 
|VB| « 31.1 G/cm. The BEC is adiabatically loaded from the 
trap into a cubic 3D optical lattice generated by three mutu- 
ally perpendicular retro-reflected laser beams at a wavelength 
of A] = 1064.5 nm to induce the phase transition to a 3D Mott 



insulator The final value for the lattice depth Vi in each direc- 
tion (i=x,y^z) is typically 20Er. For some of the experiments 
the depth is set to a higher value. Here, Er = 1.325 kHz EOl 
denotes the photon recoil energy for Cs atoms at X\. We take 
care to prepare the system with a clean singly occupied Mott 
shell with less than 4% of residual double occupancy. When 
the dipole trap is extinguished the residual harmonic confine- 
ment in the direction of gravity as a result of the transversal 
profile of the lattice beams is v, = 1 1.9(0.2) Hz. A tilt up to 
£■ = 1.7 kHz along the vertical direction can be introduced by 
lowering the strength of the levitation field. A broad Feshbach 
resonance allows us to precisely tune the atomic scattering 
length Os and thus to set U IIITI 12211 . 

We drive the phase transition from the PM to the AFM state 
in parallel for an ensemble of about 2000 ID systems ("tubes" 
or "Ising chains") with an average length of 40 sites Il23ll . We 
start with a pre-tilt of £ w Q.IU and then lower V^ to IOEr 
within 10 ms to allow for particle motion along the vertical 
direction. Effectively, our 3D one-atom-per-site Mott insula- 
tor has now turned into an ensemble of ID tubes with near 
unity occupation as horizontal tunnel coupling between the 
tubes can be neglected on the timescale of the experiment. We 
ramp through the PM-to-AFM transition with a ramp speed of 
E = 17.3 Hz/ms up to a final value Ei and project the ensem- 
ble of ID many-body systems onto number states by quickly 
ramping Vj up to its initial value and instantly reversing the 
tilt. We detect double occupancy through Feshbach molecule 
formation and detection |23 -25] with an overall doublon de- 
tection efficiency of 80(3)%. Fig. fTlplots the number of dou- 
blons A^d as a function of Ef. The transition from an ensem- 
ble of singly occupied sites to an ensemble with predominant 
double occupation can be seen clearly. The transition point 
is in very good agreement with the expectation from the BH 
model, with U = 0.963(20) kHz and E^ = 1.010(20) kHz. 
When ramping fully through the transition we create Nii — 1> 
to 4 X 10* doublons, corresponding to ensembles with 79(2)% 
double occupancy. We believe that defects in the initial one- 
atom Mott shell due to finite temperature and the superfluid 
outer ring are the main limitations to a higher conversion ef- 
ficiency. We estimate the average length of AFM chains to 
« 13 lattice sites |23|. Note that this state with exceptionally 
high double occupancy is ideal for creating ultracold samples 
of ground-state molecules |25|. 

An unambiguous characteristic of a phase transition is its 
reversibility. Fig. l2la) shows the result when the linear ramp 
(from to 40 ms) is reversed (from 40 to 80 ms). We find that 
a large fraction of the doublons, more than 85%, is returned to 
the one-atom-per-site Mott phase. 

Adiabatic driving of the quantum phase transition from the 
PM into the AFM state is permitted by the finite tunneling rate 
J between the sites through quantum fluctuations. In Fig.l2|b) 
we drive the phase transition at various values for the lattice 
depth Vj in order to highlight the role of tunneling. To assure 
adiabaticity, we increase the ramp duration for larger V^. Evi- 
dently, a reduced tunnel coupling for larger values of V^ gives 
rise to a narrowing of the transition. Note that the transition 
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FIG. 3: (color online) Dynamic response for a quench to the reso- 
nance point E ~ U with U = 1.019(20) kHz. (a) Number of dou- 
blons A'jj as a function of hold time t^ for £ = 1 .038 (squares), 0.973 
(triangles), and 0.865 kHz (circles). The solid lines are exponen- 
tially damped sinusoids fit to the data, (b) Extracted frequency fg as 
a function of E. The dashed line is a quadratic fit to the data in the 
vicinity of the minimum, (c) Extracted damping rate y as a function 
of E. The dashed line is a Lorentzian fit including a finite offset, 
(d) Extracted number of doublons A?™" at the first oscillation maxi- 
mum as a function of E. The dashed line is a Lorentzian fit. For all 
measurements V, = IOEr, Vx.y = 20 Er, and Os = 238(5) ag. 



point is shifted to a larger tilt as V, is increased due to the in- 
crease in U. We fit each data set with an error function of the 
form oc ed{{E — Eo)/w) and plot the extracted half-width w 
as a function of V, in Fig. l2tc). We clearly find the width of 
the phase transition to scale with the tunnel matrix element. 
Interestingly, w fits very well to the value 47. We suspect this 
value to be closely related to many-body enhanced tunneling 
along the chain. This is discussed next by studying quench 
dynamics in the Ising chains. 

Slow ramping through the critical point allows the system 
to follow adiabatically the many-body ground state. Non- 
equilibrium dynamics, however, arises when the initial Mott 
state is quenched by a rapidly applied tilt to the vicinity of the 
transition. In our experiment we initiate non-equilibrium dy- 
namics by first tilting the lattice and then quickly lowering V^ 
to IOEr. We have checked that the subsequent dynamic re- 
sponse does not depend on the precise way how we approach 
the transition, in particular on the order of lattice and mag- 
netic field ramps. We let the system evolve for a hold time t^ 
and finally record the number of doublons Nd as before. For 
the subsequent measurements U = 1.019(20) kHz. Fig.lsja) 
shows the result for different values of E close to the transi- 
tion. The ensemble of ID chains exhibits large amplitude os- 
cillations for A^d- The amplitude is maximal for 5 ^E — U ~Q 
and significantly lower for 1 5 1 > 0. For 5 w 0, which we call 
the resonant case, the first maximum is reached at fh w 3 ms. 



The oscillations damp out to an equilibrium value over a time 
scale of about 30 ms. We fit the data using exponentially 
damped sinusoids. From these we obtain the frequency fs, 
the peak value of the first oscillation A^™^", and the decay rate 
Y- Away from resonance, Af^"'"' is lower, however fg is higher 
and, interestingly, y is lower. This can be seen from Fig. [3jb), 
(c), and (d), where we plot fs, 7, and A^™"" as a function of 
E. We fit fg by a quadratic function, as will be justified be- 
low, and obtain the position of the minimum E^. The data for 
^max ^jj^ y ^g fitted by Lorentzians. The damping is con- 
sistent with collapse of oscillations in coherent dynamics as 
observed in our simulations of the constrained Ising model, 
see also Refs. 113] [T6i . We find that the positions of the max- 
ima for A/^"^" and 7 at 1035(2) and 1030(3) Hz, respectively, 
agree with the minimum for fg at E^ = 1028(1) Hz. Within 
our errors bars |26|, all values are compatible with the value 
for U. We note that A^^"" is only about 5% lower than the 
value for A^d we obtained from the adiabatic sweep as shown 
in Fig. [1] 

The quadratic fit to the data for fg is justified, in the sim- 
plest approximation, by a two-state Rabi model for the case 
of a single double-well potential 1(1411 . At the point E k,U, 
the two Fock states 1 1 1) and |20), describing one atom in each 
site and two atoms in the left site, respectively, are resonantly 
coupled by 7, leading to Rabi oscillations between the two 
states with frequency \/27, where the factor ^/2 ensues from 
bosonic enhancement. The probability to find the system in 
the dipole state |20) thus oscillates in the case of finite de- 
tuning from resonance with the generalized Rabi frequency 
fg = \/87^ + 5^ lfT4l . which is quadratic for sufficiently small 
5. 

Before we continue this discussion we make use of our ca- 
pability to tune U by means of a Feshbach resonance. We 
repeat the measurements shown in Fig. [3] for different values 
of fls to which we tune after producing the 3D one-atom Mott- 
insulator state and before initializing the quench. We deter- 
mine the resonance position E^ and the on-resonance oscilla- 
tion frequency /o — fs=Q- The result is shown in Fig.|4|a) and 
(b). We find good agreement for Ei- with the prediction of the 
BH model for repulsive interactions. For attractive interac- 
tions there is a significant shift to lower absolute values of E^ 
as corrections taking into account higher bands and the prop- 
erly regularized pseudopotential would have to be included 
II22I . We attribute the considerable dependence of the data for 
/o on flfs to an effectively altered tunnel barrier as U and hence 
El- are changed 1 271 1281 . 

We now quantitatively compare our findings for /o with the 
result from the simplified double-well system. Firstly, to re- 
move the effects of the modified tunneling rate due to the tilt 
and interactions, we interpolate our data to find /o for Og — 
with the linear fit as shown in Fig. Wb). The resulting val- 
ues are shown in Fig. |4|c) for data sets taken at three dif- 
ferent values of V^, and are compared to the prediction for 
the double well on resonance, /o = 2\^J (dotted line). Ev- 
idently, the tunneling in the ID chains is significantly faster 
when compared to the double-well system. This arises due 
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In conclusion, we have investigated quench dynamics by 
tuning a tih suddenly onto the transition from PM to AFM 
Ising order in a large ensemble of ID atomic Mott-insulator 
chains. We observe few-body shifts in the tunneling param- 
eter arising from interactions, and many-body effects on the 
oscillations arising from the effective constraint in the model. 
This study of quench dynamics opens the opportunity to ex- 
plore many aspects of the dynamics in these systems that up 
to now were only addressed theoretically. This includes scal- 
ing relations for smaller quenches across the phase transition 
llT5ll29l[30l . as well as potentially quenches at other resonance 
points (e.g., near E = t//2). 
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tions to numerical code development. We gratefully acknowl- 
edge funding by the European Research Council (ERC) un- 
der Project No. 278417, and support in Pittsburgh from NSF 
Grant PHY- 1148957. 



FIG. 4: (color online) Dependence of the resonance position E^ (a) 
and of the on-resonance oscillation frequency /o (b) on a^ as deduced 
from the dynamical response of the system. The solid line in (a) is 
the linear prediction of the BH model. In (b) the solid line is a linear 
fit to the data. For these measurements, Vj = IOEr and V^.y = 20Er. 
(c) On-resonance oscillation frequency /o, interpolated to a^ = 0, 
as a function of V, (circles). The dashed and dotted lines show the 
expected result for an Ising chain of three sites and the double well 
in the non-interacting limit, respectively, (d) On-resonance quench 
dynamics for the number of doublons in the ID Ising chain (circles) 
and the number of singly occupied sites in nearly isolated double 
wells (triangles) for U = 1.060(20) kHz. 



to the nearest-neighbor constraint present in the ID chains, 
which prohibits having two adjacent dipoles and which leads 
to an increase of the oscillation frequency even for three sites 
by a further factor of V2, to 47 (dashed line) lfT4l . Through 
numerical simulations with up to 60 lattice sites we have ver- 
ified that the oscillation frequency shows only a slight further 
increase for longer chains beyond 3 sites. The weak harmonic 
trapping potential causes an additional small increase of the 
computed oscillation frequency. However, both effects are 
presently within our experimental uncertainty. 

Finally, we compare the on-resonance oscillatory dynamics 
of the ensemble of ID Ising chains with that of an ensem- 
ble of nearly isolated effective double-well systems to check 
the role of the constraint. The double wells are engineered 
by highly non-adiabatic loading of the 3D lattice, thereby cre- 
ating a comparatively dilute sample of doublons. We clean 
out the singly occupied lattice sites by temporarily hiding the 
doublons in a molecular state and removing singular atoms by 
resonant light. We subsequently initiate the quench dynamics 
as before and detect the number of singly occupied sites Il23l 
as a function of fh as shown in Fig. |4|d). The oscillation fre- 
quencies in the two systems, determined to /o = 128(2) Hz 
and /o — 93(2) Hz, respectively, differ by a factor 1.38(3), in 
good agreement with the many-body effect discussed above. 
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FIG. 5: (color online) (a) Typical lattice modulation spectra showing 
tunneling resonances at modulation frequencies £ ± f/ as detected 
by an increase in the number of doublons A'd- We have set as = 
50 ao and E = 1.73 kHz (squares), a^ = 100 ao and E = 1.73 kHz 
(triangles), and as = 50 ao and E = 0.69 kHz (circles). The solid lines 
are Lorentzian fits to determine the center position of the resonances, 
(b) Tilt per lattice site £ as a function of the fraction of the magnetic 
levitation field. The dashed line is a linear fit to the data. The data 
is taken at a^ = 50 ao- (c) On-site interaction energy (/ as a function 
of as. The prediction from the BH model is depicted by the dashed 
line. In (b) and (c) all error bars are smaller than the data points. 



SUPPLEMENTARY MATERIAL: MANY-BODY QUANTUM 
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FIG. 6: (color online) (a) Calculated number of atoms as a function 
of tube length using the Gutzwiller mean field method, (b) Calcu- 
lated number of atoms as a function of chain length after introducing 
5 % randomly distributed defects. The resulting average chain length 
is calculated to ~ 13 lattice sites. 



in Fig. |5jc). The data are in very good agreement with the 
prediction from the BH model for a^ > 0, whereas for as < 0, 
U is considerably shifted to smaller absolute values due to 
multi-body interaction effects II221I . 



Estimating the average chain length 



Calibration of E and U 

The tilt E is set by reducing the vertical magnetic field gra- 
dient to a certain fraction. We calibrate E by lattice modula- 
tion spectroscopy. For the calibration measurements we set E 
to a value far detuned from the resonance point Ei- such that 
tunneling remains suppressed by the on-site interaction. Mod- 
ulation of Vj with a frequency v± ~ E±U bridges the energy 
gap and results in driven tunneling along and against the ap- 
plied potential gradient onto neighboring sites 1 1 1. 

We modulate the lattice sinusoidally for 100 ms with an 
amplitude of typically 5% around its mean value V^ = IOEr. 
Three typical excitation spectra in the number of doublons 
created are shown in Fig. l5ja). In each spectrum we identify 
two strong and narrow tunneling resonances at v± = E ±U. 
The resonance positions provide E and U. Keeping as con- 
stant they allow for precise calibration of £ as a function of 
the magnetic levitation field, as shown in Fig.lSjb). Analogous 
calibration of L'^ as a function of as for a fixed tilt is depicted 



To get a rough estimate for the average chain length for our 
ensemble of tubes we first numerically determine the ground 
state of the 3D Mott insulator using the Gutzwiller mean 
field method [2J, feeding in our typical values for the parti- 
cle number, the trapping strength, and the interaction strength. 
Fig.lSta) shows the calculated atom number distribution. The 
average atom number per tube is w 40. Defects are introduced 
by removing atoms from lattice sites at random positions. The 
remaining uninterrupted chains with a certain length can then 
be counted. Fig.l6|b) shows the resulting distribution when in- 
troducing 5 % defects. The average chain length decreases to 
PS 13 lattice sites, but evidently a reasonable number of longer 
chains survives. 



Detection of single and double occupancy 

We detect double occupancy by associating atom pairs on 
a lattice site to weakly bound molecules crossing a narrow 



g-wave Feshbach resonance with a pole at 19.8 G ["Sl. The 
system is cleaned from remaining singly occupied sites by 
combining rapid adiabatic passage from \F — 3,mf ~ 3) to 
\F =4, nif —4) using a microwave field with a resonant light 
pulse f4l. After dissociating the molecules again we detect 
the number of atoms with standard absorption imaging. Sin- 
gle occupancy is measured by associating doubly occupied 
sites into Feshbach molecules, thereby hiding these during the 
imaging process. 

Our measurement noise is mainly affected by atom num- 
ber fluctuations in the BEC. When detecting double occu- 
pancy, additional noise arises from small shot-to-shot varia- 
tions in the molecule association, cleaning, and dissociation 
efficiency. Especially for experiments with comparatively low 
absolute atom numbers, detection of single occupancy is fa- 
vored as noise from the cleaning and dissociation process is 
absent. 
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